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Abstract. Elastic-scattering phase shifts for four-nucleon systems are studied 
in an ab-initio type cluster model in order to clarify the role of the tensor force 
and to investigate cluster distortions in low energy d+d and t+p scattering. In 
the present method, the description of the cluster wave function is extended 
from a simple (0s) harmonic-oscillator shell model to a few-body model with 
a realistic interaction, in which the wave function of the subsystems are de- 
termined with the Stochastic Variational Method. In order to calculate the 
matrix elements of the four-body system, we have developed a Triple Global 
Vector Representation method for the correlated Gaussian basis functions. 
To compare effects of the cluster distortion with realistic and effective in- 
teractions, we employ the AV8' potential as a realistic interaction and the 
Minnesota potential as an effective interaction. Especially for 1 5o, the calcu- 
lated phase shifts show that the t+p and h+n channels are strongly coupled 
to the d+d channel for the case of the realistic interaction. On the contrary, 
the coupling of these channels plays a relatively minor role for the case of the 
effective interaction. This difference between both potentials originates from 
the tensor term in the realistic interaction. Furthermore, the tensor interac- 
tion makes the energy splitting of the negative parity states of 4 He consistent 
with experiments. No such splitting is however reproduced with the effective 
interaction. 

1 Introduction 



The microscopic cluster model is one of the successful models to study the struc- 
ture and reactions of light nuclei pQ. In the conventional cluster model, one 
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assumes that the nucleus is composed of several simple clusters with A < 4 
which are described by (Os) harmonic-oscillator shell model functions, and use 
an effective N-N interaction which is appropriate for such a model space. How- 
ever, it is well known that the ground states of the typical clusters d, t, h and 4 He 
have non-negligible admixtures of D-wave component due to the tensor inter- 
action. Since the conventional cluster model does not directly treat the D-wave 
component, the strong attraction of the nucleon-nucleon interaction due to the 
tensor term is assumed to be renormalized into the central term of the effective 
interaction. 

Recently, ab-initio structure calculations [2] have been successfully developed: 
Stochastic Variational Method (SVM) O HI EJ [6] , Global Vector Representation 
method (GVR) [7J [8], Green's function Monte Carlo method [9], no core shell 
model [10], correlated hyperspherical harmonics method [UJ, unitary correlation 
operator method [12j . and so on. Although the application of ab-initio reaction 
calculations with a realistic interaction are restricted so much in comparison 
with structure calculations, it has been intensively applied to the four-nucleon 
systems t+n and h+p O 13 [HI [HI [HI [[9] . Especially d+d scattering states, 
which couple to t+p and h+n channels, have attracted much attention, because 
the d+d radiative capture is one of the mechanisms making 4 He through electro- 
magnetic transitions [201 [21] and also have posed intriguing puzzles for analyzing 
powers |22[ [23| [24"1 [25| [26] , which are motivated by the famous A y problem in 
the three-nucleon system. 

Furthermore, the d+d elastic-scattering phase shifts are interesting because 
the astrophysical S-factor of the (f((f,7) 4 He reaction is not explained by any 
calculation using an effective interaction that contains no tensor term, and is 
expected to be contributed by the D-wave components of the clusters through 
E2 transitions [271 128] . 

Also, thanks to recent developments of the microscopic cluster model, the 
simple model using the (0s) harmonic-oscillator wave function with an effective 
interaction is not mandatory any more, at least, in light nuclei. We can use a kind 
of ab-initio cluster model which employs more realistic cluster wave functions 
with realistic interactions. Therefore, it is interesting to see the difference between 
the ab-initio reaction calculations with a tensor term and the conventional cluster 
model calculations without a tensor term in few-body systems. 

The microscopic i?-matrix method (MRM) with a cluster model (GCM or 
RGM) has been applied to studies of many nuclei [291 E01 EU [32] . It is now used 
in ab initio descriptions of collisions [16]. We have also applied the MRM to the 
h+p scattering problem with more realistic cluster wave functions by using a 
realistic interaction [13]. The Gaussian basis functions for the expansion of the 
cluster wave functions are chosen by a technique of the SVM pj]. In the MRM, 
as will be shown later, the relative wave function between clusters (a and b) is 
connected to the boundary condition at a channel radius. The problem is how 
to calculate the matrix elements. In this paper, we develop a method called the 
Triple Global Vector Representation method (TGVR), by which we calculate 
the matrix elements in a unified way. Although we restricted ourselves to four 
nucleon systems in the present paper, the formulation of the TGVR itself can 
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be applied to more than four-body systems as in the previous studies of the 
Global Vector Representation methods (GVR) [8|. Furthermore, for scattering 
problems, the TGVR can deal with more complicated systems than the double 
(or single) global vector which was given in the previous papers [8], because 
we need three representative orbital angular momenta, the total internal orbital 
momenta of both clusters and the orbital momentum of their relative motion, 
in order to reasonably describe the scattering states. In other words, the first 
global vector represents the angular momentum of cluster a, the second global 
vector represents the angular momentum of cluster b, and the third global vector 
represents the relative angular momentum between the clusters. 

In this paper, we will investigate the effect of the distortion of clusters on the 
d+d elastic-scattering by comparing the phase shifts calculated with a realistic 
and an effective interaction. In section 2, we explain the MRM in brief. In section 
3, the correlated Gaussian (CG) method with the TGVR, which has newly been 
developed for the present analysis, will be presented. In section 4, we will explain 
how to calculate the matrix elements with TGVR basis functions. The typical 
matrix elements are also given in the appendix. In section 5, we will present and 
discuss the calculated scattering phase shifts in detail. Finally, summary and 
conclusions are given in section 6. 



2 Microscopic i?-matrix method 

In the present study we calculate d+d and t+p (and h+n) elastic scattering phase 
shifts with the microscopic Z?-matrix method. Though the method is well docu- 
mented in e.g. Refs. [29j EH E2]) we briefly explain it below in order to present 
definitions and equations needed in the subsequent sections. Since our interest 
is on low-energy scattering, we consider only two-body channels. A channel a 
is specified by the two nuclei (clusters) a, 6, their angular momenta, I a ,h, the 
channel spin / that is a resultant of the coupling of I a and lb, and the orbital 
angular momentum £ for the relative motion of a and b. The wave function of 
channel a with the total angular momentum J, its projection M, and the parity 
7r takes the form 



I 



*?X xM , (2-1) 



JM 



where <&" a and & b Ib are respectively antisymmetrized intrinsic wave functions of 
a and 6, and A is an operator that antisymmetrizes between the clusters. The 
square bracket [I a Ib]i denotes the angular momentum coupling. The coordinate 
p a in the relative motion function Xa is the relative distance vector of the clusters. 
The channel spin / and the relative angular momentum £ in a are coupled to give 
the total angular momentum J. The relative-motion functions Xa also depend 
on J and ir. For simplicity, this dependence is not displayed explicitly in the 
notation for Xa as well as for some other quantities below. 

The configuration space is divided into two regions, internal and external, by 
the channel radius a. In the internal region (p a < a), the total wave function 



may be expressed in terms of a combination of various VP 
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_(_(_l) A a+/a+4-/+^ 



JM 



JM 



(2.3) 



where j4 a is the number of nucleons in cluster a, 6 a b is unity if a and & are identical 
clusters and zero otherwise, and 5j a j h is unity if the clusters are in identical states 
and zero otherwise. In the second line of Eq. (|2.2p . the relative motion functions 
of Eq. (|2.1h are expanded in terms of some basis functions as 



Xam(pa) ^ ^ fanX^an 

(p a )Y im (p a ). (2.4) 

n 

In what follows we take 

U a n(pa) = Pa SXp(- - \ n p a ) (2.5) 

with a suitable set of A n s. 

In the external region (p a > a), the total wave function takes the form 



JMtt 
ext 



^ g a ( Pa )<P J a M \ 



a 



(2.6) 



Note that the antisymmetrization between the clusters is dropped in the external 
region under the condition that the channel radius a is large enough. The function 
9a(Pa) of Eq. (|2,6p is a solution of the equation 



J^(_f_ 2 d £(£+!) 

2 p a \dp a p a dp a p 2 a 



+ 



Z a Z h e 7 

Pa 



g a (p a ) = E a g a (p a ), (2.7) 



where p a is the reduced mass for the relative motion in channel a, Z a e and Z^e 
are the charges of a and b, and E a = E — E a — is the energy for the relative 
motion, where E is the total energy, and E a and Ef, are the internal energies 
for the clusters a and b, respectively. For the scattering initiated through the 
channel ckq, the asymptotic form of g a for the open channel a {E a > 0) is 



9a(Pa) = v a 

Pa [Ia(kaPa)fiaao S aa O a (k a p a )] 



(2.8) 



where k a = y/2p, a \E a \/h, v a = hk a /p a and is an element of the 5-matrix 
(or collision matrix) to be determined. Here I a (kaPa) and O a (k a p a ) are the 
incoming and outgoing waves defined by 



la(kaPa) = O a (k a p a )* = Gi(r] a ,k a p a ) - iFe(r) a ,k a p a ), 



(2.9) 
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with the regular and irregular Coulomb functions F% and Gg. The Sommerfeld 
parameter r\ a is p a Z a Zi,e 2 /h 2 k a . For a closed channel a (E a < 0), the asymptotic 
form of g a is given by the Whittaker function 

9a{Pa) OC P^W^^^kaPa). (2.10) 

The matrix elements S^ are determined by solving a Schrodinger equation 
with a microscopic Hamiltonian H involving the A a + nucleons, 

(H + C - E)*ttf* = (2.11) 

with the Bloch operator £ 

where the channel radius a is chosen to be the same for all channels, and the b a 
are arbitrary constants. Here, we choose b a = for the open channels and b a = 
2k a aW'_ r] £ +1 / 2 (2^a a )/W / -j7 a /+i/2(2^a«) for the closed channels. The results do 
not depend on the choices for b a but these values simplify the calculations. 
Notice that the projector on \(f>i M7T ) in Eq. (|2.12p is not essential in a microscopic 
calculation and can be dropped since the various channels are orthogonal at the 
channel radius. 

The Bloch operator ensures that the logarithmic derivative of the wave func- 
tion is continuous at the channel radius. In addition, ^i^ In must be equal to 
^cxt n & t Pa = a,. Projecting the Schrodinger equation on a basis state, one 
obtains 

E C -'n',an fan = (^\C\^) (2.13) 



an 



with 



and 



CWan = {* J Jt?\H + C- E\A$ J J? n )- mt , (2.14) 



$ J J?* = u an (p a )4> J a M \ (2.15) 

Here (|0|)mt indicates that the integration with respect to p a is to be carried 
out in the internal region. Actually (|C|)mt is obtained by calculating the ma- 
trix element (\0\) in the entire space and subtracting the corresponding external 
matrix element (|0|) cx t that is easily obtained because no intercluster antisym- 
metrization is needed. The i?-matrix and Z-matrix are defined by 



( k a / \ 2 % > 1\ / \ / \ 

= -5- ( 7— I 2_^Ua'n'{a){C ) a 'n',anUan{a), (2.16) 

^ \pa' pa^a / , 
n'n 

Za'a = Ia{k a a)5 a i a - TZ a i a k a al' a (k a a) . ( 2 -17) 
The S'-matrix is finally obtained as 

S Jw = {Z*)~ l Z. (2.18) 
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Table 1. Channel spins ( 2I+1 ij) of physical d+d, t+p, and h+n channels for J < 2 and I < 2. 
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3 Po 


'Pi 
3 Pi 


3 P2 



In this paper we focus on the elastic phase shifts 8^ that are defined by the 
diagonal elements of the S-matrix, 

Stl = V J Je™l\ (2.19) 

We study four-nucleon scattering involving the d+d, t+p and h+n channels 
in the energy region around and below the d+d threshold. In Table [1] we list all 
possible labels 21+1 1 j of physical channels for J n = 0^, 1^, and 2^, assuming £ < 
2. Here "physical" means that the channels involve the cluster bound states that 
appear in the external region as well. Non-physical channels involving excited 
pseudo states will also be included in most calculations. Note that the d+d 
channel must satisfy the condition of I + £ even (see Eq. (|2.3p ). The channel spin 
/ = or 2 can couple with only even £, but 1 = 1 with only odd £. It is noted 
that the relative motion for the d+d scattering can have I = only when J 11 is 
equal to + and 2 + . 

Because one of our purposes in this investigation is to understand the role of 
the tensor force played in the four-nucleon dynamics, we want to compare the 
phase shifts obtained with two Hamiltonians that differ in the type of NN inter- 
actions. One is a realistic interaction called the AV8' potential [33] that includes 
central, tensor and spin-orbit components. We also add an effective three-nucleon 
force (TNF) in order to reproduce reasonably the binding energies of t, h and 
4 He |34j . which makes reasonable thresholds. In the present calculation, the TNF 
is included in all calculations for AV8'. Another is an effective central interaction 
called the Minnesota (MN) potential [35], which reproduce reasonably the bind- 
ing energies of t, h and 4 He, though it has central terms alone (with an exchange 
parameter u = 1). The Coulomb potential is included for both potentials. 

The intrinsic wave function of cluster k (k = a, b) is described with a 
combination of basis functions with different and Sk values 

*u fc =J2 a Kr c) 4 s r n) ] IhM 4:°m? ) . ( 2 - 2 °) 



where ip^^ , ^Sk^ anc ^ ^T^M^ denote the space, spin and isospin parts of the 



cluster wave function. In the case of the AV8' potential, the t (or h) wave func- 
tion is approximated with thirty Gaussian basis functions that include L}, < 2, 
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Table 2. Energies E, rms radii i? rms and D-state probabilities Pd of the clusters that appear 
in four-nucleon scattering and 4 He with the AV8' (with TNF) and MN potentials. Nk is the 
number of basis functions used to approximate the wave function of cluster k. The values in 
the last three columns for three- and four-body systems are taken from Ref. [34] for AV8' and 
Ref. [7] for MN. 



potential 
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present 




literature 








N k 


E 


rjrms 


Pd 


E 




Pd 








(MeV) 


(fm) 


(%) 


(MeV) 
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(%) 






8 


-2.18 


1.79 


5.9 


-2.24 


1.96 


5.8 


AV8' 




30 


-8.22 


1.69 


8.4 


-8.41 






(with TNF) 


4 He(0+) 


30 


-7.55 


1.71 


8.3 


-7.74 








(2370) 


-27.99 


1.46 


13.8 


-28.44 




14.1 




*(D 
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-2.10 


1.63 





-2.20 


1.95 





MN 


15 


-8.38 


1.70 





-8.38 


1.71 







4 He(0+) 


15 


-7.70 


1.72 





-7.71 


1.74 







(1140) 


-29.94 


1.41 





-29.94 


1.41 






and Sk = 2 an d f • The deuteron wave function is also approximated with Gaus- 
sian basis functions, four terms both in the S- and D- waves, respectively. The 
falloff parameters of the Gaussian functions are selected using the SVM [5] and 
the expansion coefficients are determined by diagonalizing the intrinsic cluster 
Hamiltonian. A similar procedure is applied to the case of the MN potential. 

The calculated energies E, root-mean-square (rms) radii R rms and D state 
probabilities Pd are given in the fourth to sixth columns in Table El We use the 
truncated basis in order to obtain the phase shifts in reasonable computer times, 
they slightly deviate from more elaborate calculations, which are given in the last 
three columns. Fortunately, except for the small shift of the threshold energy, 
the phase shifts are not very sensitive to the details of the cluster wave functions 
because they are determined by the change of the relative motion function of the 
clusters. The Nk values in parenthesis for 4 He are the number of J T = + basis 
functions in the major multi-channel calculation. The energy of 4 He calculated 
in Table 2 with the multi-channel calculation is thus not optimized but found 
to be very close to that of the more extensive calculation. It is noted that the 
calculated i? rms value for the deuteron is smaller than in other calculations. This 
is due to the restricted choice of the length parameters of the basis functions, 
which permits us to use a relatively small channel radius of a ~ 15 fm. We have 
checked that the phase shifts for a d + d single channel calculation do not change 
even when more extended deuteron wave functions are employed. 
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3 Correlated Gaussian function with triple global vectors 



As explained in the previous section, the calculation of the S'-matrix reduces to 
that of the Hamiltonian and overlap matrix elements with the functions defined 
by (|2.ip and (|2.20p . and it is conveniently performed by transforming that wave 
function into an LS'-coupled form, 



A 



^ (space) ^ (space) 



Xa(pa) 



L 



./.(spin) / (spin) 



(3.1) 



JM 



The transformation can be done as 

^[[[vAr cc v? in) ]/j< pacc vi pin) ] 



XaiPc 



JM 



, (isospin) , (isospin) 
^T a M Ta Vt o M t 



L a bLS 

xA 



L a S a I a 
Lb Sb lb 
L a b S I 



(space) 

L b (L ab )£L 



i/, (spin) 7// spin) 



tfak+J-i-LutSLrtJl-JL) 

(isospin) , (isospin) 



JM V T a M Ta V Tb M Tb 



with 



, (space) 
^L a L b (L ab )eL 



^(space)^(spacc) 



Xa(p a ) 



(3.2) 



(3.3) 



where U and [ ] are Racah and 9j coefficients in unitary form [5]. 

The evaluation of the matrix element can be done in the spatial, spin, and 
isospin parts separately. The spin and isospin parts are obtained straightfor- 
wardly. In the following we concentrate on the spatial matrix element. The spa- 
tial part (|3.3p of the total wave function is given as a product of the cluster 
intrinsic parts and their relative motion part. The coordinates used to describe 
the 2N+2N channel are depicted in Fig.[I]|a) with p a = X3, whereas the coordi- 
nates suitable for the t+p and h+n channels are shown in Fig. HJb) with p a = x' 3 . 
These coordinate sets are often called H-type and K-type. Therefore the calcula- 
tion of the spatial matrix element requires a coordinate transformation involving 
the angular momenta L a ,Lb,i and L a b- Moreover the permutation operator in 
A causes a complicated coordinate transformation. All these complexities are 
treated elegantly by introducing a correlated Gaussian [3(3 [H [5] , provided each 
part of tb^T V r CC /\ \nr is given in terms of (a combination of) Gaussian functions 

r L a L b (L ab )ZL . 

as in the present case. In what follows we will demonstrate how it is performed. 
Because the formulation with the correlated Gaussian is not restricted to four 
nucleons but can be applied to a system including more particles, the number of 
nucleons is assumed to be N in this and next sections as well as in Appendices 
B and C unless otherwise mentioned. 

The relative and center of mass coordinates of the N nucleons, Xi (i = 
1,...,N), and the single-particle coordinates, r j (i = 1,...,N), are mutually 
related by a linear transformation matrix U and its inverse U~ l as follows: 



x \ 



N 

.7=1 



A' 



(3.4) 







3 1 



(a) H-type 



4 2 




(b) K-type 



Figure 1. Relative coordinates for the four-body system 



We use a matrix notation as much as possible in order to simplify formulas and 
expressions. Let x denote an (N — l)-dimensional column vector comprising all 
Xi but the center of mass coordinate xn- Its transpose is a row vector and it is 
expressed as 

x = (x 1 ,x 2 , ...,£Civ-i). (3.5) 

The choice for x is not unique but a set of Jacobi coordinates is conveniently 
employed. For the four-body system, the Jacobi set is identical to the K-type 
coordinate, and the corresponding matrix U is given by 



U K 



( 1 

l 

2 
1 

? 
1 



-1 

\ -1 

1 1 

? I 

4 4 





-1 

I J 



(3.6) 



The transformation matrix for the H-type coordinate reads 



( 1 



i 

2 

VI 





1 

1 

' 2 
I 

1 



\ 
-1 
1 

"2 



(3.7) 



The K-type coordinate is obtained directly from the H-type coordinate by a 
transformation matrix UkU 



h 



UkUh 1 



( 1 



V o 






1 / 



Ukh 
1 



(3.8) 



where Ukh is a 3x3 sub-matrix of UkU J] . 

Each coordinate set emphasizes particular correlations among the nucleons. 
As mentioned above, the H-coordinate is natural to describe the d+d channel, 
whereas the K-coordinate is suited for a description of the 3A^+A^ partition. 
It is of crucial importance to include both types of motion in order to fully 
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describe the four- nucleoli dynamics [2]. In order to develop a unified method 
that can incorporate both types of coordinates on an equal footing, we extend 
the explicitly correlated Gaussian function (37J E] to include triple global vectors 

f l 1 l 2 (L 12 )L 3 lm(ui,U2,u 3 , A, x) 

= exp (-^xAxj [[y Ll (^x)y L2 (l^x)} Ll2 y L3 (^x)} LM , (3.9) 



where 



yLiMi(uix) = \uix\ Ll Y LiMl (uix) (3.10) 



is a solid spherical harmonics and its argument, U{X, what we call a global vector, 
is a vector defined through an (N — l)-dimensional column vector itj and x as 

N-l 

UiX = ^2(v,i)jXj, (3.11) 
i=i 

where (ui)j is the jth element of u%. In Eq. (I3.9P A is an (N — 1) x (N — 1) 
real and symmetric matrix, and it must be positive-definite for the function F 
to have a finite norm, but otherwise may be arbitrary. Non-diagonal elements of 
A can be nonzero. 

The matrix A and the vectors u\ , u 2 , u 3 are parameters to characterize the 
"shape" of the correlated Gaussian function. The Gaussian function including A 
describes a spherical motion of the system, while the global vectors are respon- 
sible for a rotational motion. The spatial function (j3.3j) is found to reduce to 
the general form (|3.9|) . Suppose that x stands for the H-type coordinate. Then a 
choice of ui=( 1,0,0), it2=(0,l,0) and ■U3=(0,0,l) together with a diagonal matrix 
A provides us with the basis function (|3.3[) employed to represent the configu- 
rations of the 2N+2N channel. On the other hand, the K-type basis function 
looks like 

exp LIx>a k A [[y L Axi)yu*2)W 2 y L3 {x'z)] LM , (3.12) 

where x' = (x[,X2,x' 3 ) is the K-coordinate set (see Fig. 1(b)) and Ak is a 3x3 
diagonal matrix. Noting that x' is equal to x' = Ukh x , we observe that the basis 
function f|3. 12|) is obtained from Eq. (|3.9p by a particular choice of parameters, 
that is, ui=(l,0,0), U2=(0,— ^,1) and ; U3=(0,|,|), and the matrix A is related to 
A K by 



A = (uiu 2 u 3 )A K \ u 2 \ = UkhA k Ukh- (3.13) 




Thus the form of the F-function remains unchanged under the transformation 
of relative coordinates. 

Note that A is no longer diagonal. The choice of a different set of coordinates 
ends up only choosing appropriate parameters for A, u\, u 2 , and u 3 . 
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It is also noted that the triple global vectors in Eq. (I3.9P are a minimum 
number of vectors to provide all possible spatial functions with arbitrary L and 
parity tt. A natural parity state with tt = (— 1) L can be described by only one 
global vector, that is, using e.g., L\ = L, L2 = 0, L12 = L, L 3 = [6j [381 EJ- To 
describe an unnatural parity state with tt = (— 1) L+1 except for 0~ case, we need 
at least two global vectors, say, L\ = L, L2 = 1, L12 = L, L 3 = [371 EJ- The 
simplest choice for the - state is to use three global vectors with L\ = L2 = 
L12 = L 3 = 1 [37J. In this way, the basis function (|3.9p can be versatile enough 
to describe bound states of not only four- but also more-particle systems with 
arbitrary L and tt. 

To assure the permutation symmetry of the wave function, we have to operate 
a permutation P on F. Since P induces a linear transformation of the coordi- 
nate set, a new set of the permuted coordinates, Xp, is related to the original 
coordinate set x as xp = Vx with an [N — 1) x (JV — 1) matrix V . As before, 
this permutation does not change the form of the F-function: 



The fact that the functional form of F remains unchanged under the permutation 
as well as the transformation of coordinates enables one to unify the method of 
calculating the matrix elements. This unique property is one of the most notable 
points in the present method. 

4 Calculation of matrix elements 

Calculations of matrix elements with the correlated Gaussian F are greatly fa- 
cilitated with the aid of the generating function j PU 



PF Ll L 2 (L 12 )L :i LM(ui,U2,U 3 , A, x) 

= Fl 1 l 2 (l 12 )l 3 lm(ui,U2,U3, A, X P ) 

= Fl 1 l 2 (l 12 )l 3 lm(f'u 1 ,Vu2,Vu 3 ,VAV,x). 



(3.14) 




(4.1) 



with s = (si, S2, ■ ■ ■ , Sjv-i)j where Si = Yl^=i ^j(, u j)i e j> e j 1S a 3-dimensional 
unit vector (e,- • ej = 1), and Xj is a scalar parameter. More explicitly 



iV-l 



N-l 3 



3 




(4.2) 



The correlated Gaussian F is generated as follows: 



F] 



Y 1 L 2 (L 12 )L 3 LAl{ui,U 2 ,U 3 ,A, X) 




(4.3) 
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where 

s L - <™ m 

When g(s;A,x) is expanded in powers of Ai, only the term of degree A^ 1 con- 
tributes in Eq. (|4.3p . and this term contains the Lith degree e\ because Ai and 
ei always appear simultaneously. In order for the term to contribute to the in- 
tegration over e\, these L\ vectors e\ must couple to the angular momentum 
L\ because of the orthonormality of the spherical harmonics Yl 1 m 1 (gi), that is, 
they are uniquely coupled to the maximum possible angular momentum. The 
same rule applies to A2, e 2 and A3, as well. 

We outline a method of calculating the matrix element for an operator O 

( f l 4 l 5 (l 45 )l 6 l'M'{u4, u 5 ,u G ,A', x)\0\F LlL2{Ll2)LsLM (u 1 ,u 2 , u 3 , A, x)) . (4.5) 

In what follows this matrix element is abbreviated as {F'\0\F). Using Eq. (14.3|) in 
Eq. (|4.5j) enables one to relate the matrix element to that between the generating 
functions: 

{F'\0\F) = (nff / [[^ 4 (e4m 5 (e 5 )] i45 F i6 (e 6 )]^ ; 
x [[Y Ll ( ei )Y L2 (e 2 )] Li2 Y La (e 3 )] LM 

x \Il ^ J (9(s', A', x')\0\g(s, A, x)) \^ (4.6) 



with 



s = Aiiiiei + A2U 2 e 2 + A3n3e3, s' = A4U4e4 + A 5 u 5 e5 + A 6 u 6 e 6 . (4.7) 

The calculation of the matrix element consists of three stages: (1) Evaluate the 
matrix element between the generating functions, {g(s' , A' ,x')\0\g(s, A,x)) . (2) 
Expand that matrix element in powers of Aj and keep only those terms of degree 
Li for each i. (3) Recouple the vectors ej and integrate over the angle coordinates. 
In the second stage the remaining terms should contain e^s of degree Lj as well. 
Hence any term with A?ej-ej = A? etc. can be omitted because the degree of e- L 
becomes smaller than that of Aj. 

We will explain the above procedures for the case of an overlap matrix ele- 
ment. The matrix element between the generating functions is 




with 



6 

B = A' + A, Z = S + s' = ^ ^i e i u i- 



(4.9) 
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To perform the operation in the second stage we note that 

6 



1 1 6 

-zB~ x z = - ^ pijXiXjei ■ ej (4.10) 



with 



Pij = UiB 1 u j . (4.11) 
As mentioned above, here we can drop the diagonal terms, A?e; • e^, and we get 



nSl) (9(s',A',x%(s,A,x 



A l= 



Here the non-negative integers riij must satisfy the following equations in order 
to assure the degree Lj for in the different terms, 

"12 + "-13 + "-14 + "15 + "16 = Ll, 
"12 + "-23 + "-24 + "25 + "26 = L 2 , 
"13 + "23 + "34 + "35 + "36 = L 3 , 
"14 + "-24 + "-34 + "45 + "46 = ^4 5 
"15 + "-25 + "35 + "45 + "56 = L 5 , 

"16 + "-26 + "-36 + "46 + "56 = L 6 . (4.13) 

The last step is to recouple the angular momenta arising from the various 
terms. Since we have to couple ejS to the angular momentum Lj from the terms 
of degree Li, we may replace the term (pijei ■ ej) niJ with just one piece 



1 J> B JV 3 [Yn i3 {ei)Y nij (e,-)] 00 . (4.14) 

riij 

Other pieces like [Y K (ei)Y K (ej)] 00 with k < riij do not contribute to the matrix 
element. We thus have a product of 15 terms of [Y nij (ei)Y nij (e^)] 0Q . The coupling 
of these terms is done by defining various coefficients that are all expressed in 
terms of Clebsch-Gordan, Racah, and 9j coefficients. For example, we make use 
of the formulas 

[^(ei)Y o (e 2 )] 00 [n(ei)n(e 3 )]oo [Y c (e 2 )y c (e 3 )]oo 

X(abc) [[Y a+fe (ei)Y a+c (e2)] b+c y 6+c (e 3 )]oo, (4.15) 

[Y a (e 1 )Y a {e 4 )] Q0 [Y b ( ei )Y b (e 5 )) 00 [Y c ( ei )Y c (e 6 )] 00 
-> R 3 (abc) [Y a+b+c ( ei ) [[F a (e 4 )F 6 (e 5 )] a+6 F c (e 6 )] a+6+c ]oo. (4.16) 



Here the symbol — > indicates that no other terms arising from the left hand side 
of the equation contribute to the integration over the angles e^s, so that only the 
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term on the right hand side has to be retained. Another coefficient is 

[[[Y a (e A )Y b (e b )} q Y c {e & )] Q [[Y a ,(e 4 )Y bl (e 5 )] q ,Y c ,(e 6 )] Q/ ]t 
-> W(abcqQ, a'b'c'q'Q', ££') [[Y a+a , (e 4 )Y b+b , {e b )] e Y c+c , (e 6 )]/. (4.17) 

Expressions for the coefficients, X, R 3 , W, are given in Appendix A. Performing 
the integration of the six unit vectors, e^s, as prescribed in Eq. (|4.6p leads to the 
overlap matrix element 

(F'\F) 

/ (2vr)^ \ 3/2 /« \ (-i)^+W 3 

= l^eT^J L^A 7 V2L + 1 fc£,w 

xJ2 \\h + 1 0(n 4J ;LiL 2 L3L4L 5 L 6 ,Li 2 L45L),(4.18) 

n SJ \i<i ni J / 

with 

0(iiij; L1L2L3L4L5L6, L^L^L) 

= X(ni2ni3n23)-R3(^14"-15^16)-R3("-24"-25n26)-R3(n34^35"-36)^("'45"'46"'56) 

x Z(?ii2+ni 3 L 1 -n 1 2-n 13 )Z(n 1 2+n 2 3 L 2 -ni2-n2 3 )Z(n 13 +n2 3 L 3 -n 13 -n 23 ) 

L12 £\ n 13 +n 23 

L 3 L 3 -ni 3 -n 23 n 13 +n 23 
L L 

x W(nun 15 niQ nu+n 15 Li-ni 2 -ni 3 , n 2 4^25"26 n 2 4+"-25 £2-^12-^23, ^2) 

x H^(ni4 + n 2 4 rai 5 + n 25 ni 6 + n 26 £2 4,n34«35™36 ^34+^35 L 3 -rii 3 -n 23 , Ll 3 ) 

x W^(L 4 -n45-n46 L 5 -n 45 -n 56 L e -n 4 Q- n 56 ^3 L, 

n45 + n 46 n 45 + n 56 n 46 +n 56 n 46 +n 56 0,LL 45 ), (4.19) 

where Z is the coefficient given in Eq. (lA.3j) . The summation in Eq. (I4.18D extends 
over all possible sets of n^- that satisfy Eq. (|4.13p . In most cases the values of Lj 
are limited up to 2, so that the number of terms to be evaluated is not so large 
and the calculation of the matrix element is fast. 

Expressions for the Hamiltonian matrix elements are collected in Appendix 
B. One advantage of our method is that the calculation of matrix elements can be 
done analytically. In addition we do not need to do angular momentum and par- 
ity projections because the correlated Gaussian function (|3.9p already preserves 
those quantum numbers. 

The Fourier transform of the correlated Gaussian function F is a momentum 
space function and it becomes a useful tool to calculate various matrix elements 
that depend on the momentum operators [7]. For example, the distribution of 
the relative momentum is obtained by the expectation value of 8(p i — pj — p), 
where Pj is the momentum of the jth particle. It is obviously much easier to 
calculate the distribution using the momentum space function rather than the 
coordinate space function. We show in Appendix C that the Fourier transform 
of F reduces to a linear combination of Fs in the momentum space. 



Li L 1 -n X 2-n 13 n l2 +ni 3 
^ L 2 L 2 -ni2-n 23 ni2+n 23 
£ih£3 L12 l\ ni 3 + n 23 
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5 Results 



5.1 2N+2N and 3N+N channels 

In Table[TJ we gave the physical channels, d+d, t+p, and h+n. Fig. [2] displays two- 
body decay thresholds in the d+d threshold energy region. The three physical 
channels are the main channels that describe the scattering around the three 
lowest thresholds (d+d, t+p, h+n). However, the scattering wave function VPjJjt 
in the internal region should contain all effects that may occur when all the 
nucleons come close to each other. It is thus reasonable that ^j^^ may not be 
well described in terms of the physical channels alone. Particularly the deuteron 
can be easily distorted when we use realistic potentials. 

We will show that some pseudo 2N+2N channels are indeed needed to sim- 
ulate the distortion of the deuteron. These pseudo channels, when they are in- 
cluded in the phase-shift calculation, are expected to take account of the dis- 
tortion of the clusters of the entrance channel [30] . Here "pseudo" means that 
the clusters in the pseudo channels are not physically observable but may play 
a significant role in the internal region. The wave functions of these 2N pseudo 
clusters are obtained by diagonalizing the intrinsic cluster Hamiltonian simi- 
larly to the case of the physical clusters. We take into account the following 
pseudo clusters: d*(l + ,T = 0), d*(0 + ,T = 1), d*(2+,T = 0), d*{3 + ,T = 0), 
2n*(0 + ,T = 1), and 2p*(0 + ,T = 1), where the upper suffix * indicates all the 
excited state but the ground state of d. Among the pseudo clusters, the low- 
est energy states with + that are related to virtual states would be most im- 
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2n(0+)+2p(0+) 
d(0+)+d(0+) 



d(l+)+d(0+) 



d(l+)+d(l+ 



d(l+)+d(0+ 



■ d(l+)+d(l+ 



p+p+n+n 



■ d(l+)+d(l+) 



10- 



AV8' 



h+n 
t+p 



MN 



h+n 
t+p 



Exp. 



h+n 
t+p 



Figure 2. Two-body thresholds calculated with the AV8' (left) and MN (middle) potentials. 
The solid lines are physical channels and the dashed lines are pseudo channels. We also plot 
experimental two-body thresholds for physical channels (right). The dotted line is the p+p+n+n 
threshold. 
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Table 3. 2N+2N and 3N+N channels. The Roman and Arabic numerals correspond to sets 
of channels included in the calculations. 



model 


channel 




2N+2N 


I 


d(l+)+d(l+) 








d(l+)+cf(l + ) 








d* (l+)+cf(l+) 






II 


J(0+)+J(0+) 








J(0+)+d*(0+) 








<f (0+)+d*(0 + ) 






TTT 

111 


d*(2+)+cf (1+) 








a (z^j+a (2 j 






IV 


d*(3+)+d*(l+) 


FULL 






d*(3+)+d*(2+) 








d*(3 + )+d*(3 + ) 






V 


2n(0+)+2p(0 + ) 








2n(0 + )+2p*(0+) 








2n*(0+)+2p(0+) 








2n*(0+)+2p*(0+) 




3N+N 


1 

2 


«(D+p(D 

/ l (i + )+n(i + ) 
/ 1 *(i + )+n(i + ) 



portant. We especially write them as d(0 + ), 2n(0 + ) (di-neutron) and 2p(0 + ) 
(di-proton). Although they are not bound, they are observed as resonances 
or quasi-bound states with negative scattering lengths. In fact the scattering 
lengths are a s (nn) = —16.5 fm and a s {pp) = —17.9 fm, which are comparable to 
a s (np,T = 1) = —23.7 fm. The calculated thresholds of these pseudo channels 
are also drawn in Fig. [2 

Though it is expected that the pseudo channels with low threshold ener- 
gies contribute more strongly to the scattering phase shift, we take into ac- 
count all of these 2N+2N channels that include a vanishing total isospin as 
given in Fig. [2j The total isospin of the 3iV-|-./V channel is mixed in the present 
calculation. Because the T=\ component of the scattering wave function only 
weakly couples to the d(l + ,T = 0)+d(l + ,T = 0) elastic-channel, the channel 
d(l + ,T = 0)+(f(0 + , T = 1) is not employed in the calculation. 

We also include the excited deuteron channels that comprise the d*(2 + ,T = 
0) and d*(3 + ,T = 0) clusters. The energies of these lowest thresholds are above 
10 MeV. These channels are therefore expected not to be very important, but 
that is not always the case as will be discussed in the case of the 1 So d+d phase 
shift. 

Table [3] summarizes all the channels that are used in our calculation. The 
2N+2N channels are distinguished by Roman numerals, while the 3iV-|-./V chan- 
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nels are labeled by Arabic numerals. In the following, we use an abbreviation 
"2iV+2iV" or "3iV+iV" to indicate calculations including all 2N+2N channels 
I-V or all 3iV+iV channels (1-2 in Table EJ), respectively. Here t*(\ + ) and h*(\ + ) 
are excited 3N continuum states. A "FULL" calculation indicates that all the 
channels in the table are included to set up the S'-matrix. In the case of the MN 
potential channels III and IV are not included because this potential contains 
no tensor force. 

The relative wave functions Xam are expanded with 15 basis functions. We 
checked the stability of the S'-matrix against the choice of the channel radius. 
The channel radius employed in this calculation is about 15 fm. 



5.2 Positive parity phase shifts 

Fig. [3] displays the 1 So d+d elastic-scattering phase shift obtained with the AV8' 
potential. The dash-dotted line is the phase shift calculated with channel I (Id = 
1 + ), and the dash-dot-dotted line is the phase shift with channels I and II (Id = 
l + ,0 + ). The phase shifts calculated by including further excited deuterons are 
also plotted by the dashed and dotted lines that correspond to the channels I-III 
(Id < 2 + ) and I-IV (Id < 3 + ), respectively. A naive expectation that the 1 5o d+d 
elastic-scattering phase shift might be well described in channel I (d(l + ) + d(l + ), 
d(l + ) + d* (l + ) and d*(l + ) + d*(l + )) alone completely breaks down in the case 
of the AV8' potential. 

Because the deuteron has a virtual state d with + at low excitation energy, it 
is reasonable that the inclusion of channel II gives rise to a considerable attractive 
effect of several tens of degrees on the phase shift, as shown by the dash-dot- 
dotted line of Fig. [3l However, the phase shift exhibits no converging behavior 
even when the higher spin states such as d*(2 + ) and d*(3 + ) are taken into account 
in the calculation. The additional attractions by these channels are of the same 




EcliMeV) 



Figure 3. 1 5'o d+d elastic-scattering phase shift calculated with the AV8' potential. The phase 
shifts are all obtained within the d+d channels. The set of included channels is successively 
increased from I to IV. See Table [3] for the deuteron states included in each channel. 
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order as that of channel II. One may conclude that the deuteron is strongly 
distorted even in the low energy 1 So d+d elastic scattering but more physically we 
have to realize that there exist two observed + states below the d+d threshold. 
Obviously the d+d scattering wave function is subject to the structure of those 
states in the internal region. 

The second + state of 4 He lying about 4 MeV below the d+d threshold is 
known to have a 3N+N cluster structure |39|. [34] . Thus this state together with 
the ground state of He cannot be described well in the 2N+2N model space 
alone. As seen in Table [IJ the 3N+N channel contains a 1 So component, which 
is the dominant component of the 0^ state. Since the realistic force strongly 
couples the 2N+2N channel to the 3N+N channel and the d+d scattering wave 
function has to be orthogonal to the main component of the underlying + states, 
we expect that the deuteron in the incoming d+d channel never remains in its 
ground state but has to be distorted largely due to the 3N+N channel. The 
phase shift for the channel I-IV (dotted line) shows a resonant pattern. This 
resonant state is expected to be the second + state because of the restricted 
model space within the d+d channel. 
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Figure 4. Comparison of the ground and second + state energies between calculations with 
the AV8' (left) and MN (middle) potentials and experiment (right). The model space for AV8' 
is I-IV, 3N+N and FULL and the model space for MN is I-II, 3N+N and FULL. 



Fig. [3] displays the calculated ground state energy and the second + en- 
ergy for the AV8' (left) and MN (middle) potentials. The model spaces of the 
calculations are I-IV, 3^+^ and FULL for AV8' and I-II, 3N+N and FULL 
for MN. We also plot experimental energies (right) [30] • For the AV8' potential, 
the energies of the two lowest + states do not change very much between the 
FULL and 3N+N models. But the second + state with the d+d model (chan- 
nels I-IV) is not bound with respect to the d+d threshold as expected before. 
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Figure 5. 1 S'o d+d elastic-scattering phase shift calculated with the AV8' (left) and MN 
(right) potentials. The solid line is a FULL calculation, while the dotted line is a d+d channel 
calculation. Crosses correspond to the 7?-matrix analysis of Ref. [22] . 



On the contrary, for the MN potential, the second + state with the d+d model 
(channels I-II) is bound with respect to the d+d threshold. We consider that 
this difference makes the drastic change of the d+d phase shifts, between the 
AV8' and MN potentials. It is also interesting to see that the energies of the two 
lowest + states for the MN potential are almost the same between the FULL 
and 3N+N models. 

Plotted in Fig. [5] are the 1 So d+d elastic-scattering phase shifts obtained with 
the AV8' potential (left) and the MN potential (right). The FULL calculation 
(solid line) couples all 2N+2N and 3N+N channels that are listed in Table El 
The i?-matrix analysis (crosses) |22j is reproduced well by both the AV8' and 
MN potential with the FULL calculation. Compared to the uncoupled phase shift 
(dotted line), one clearly sees that the 3N+N channel produces a very large effect 
on the d+d elastic phase shift, especially in the case of the AV8' potential. We 
also verified that a calculation excluding the channels III, IV or V from the 
FULL channel calculation gives only negligible change in the phase shift. The 
slow convergence seen in Fig. [3] is thus attributed to the neglect of the 3N+N 
channel, indicating that a proper account of the 1 Sq d+d elastic phase shift 
at low energy can be possible only when the coupled channels {d(l + )+d(l + )} 
+{d{0 + )+d(0 + )}+ {t(l/2+)+p(l/2+)} +{/i(l/2+)+n(l/2+)} are considered. 

Thus, the slow convergence in Fig. [3] suggests that the 2N+2N partition is 
not an economical way to include the effects of the 3N+N channel. In the case 
of the MN potential (right panel in Fig. [5]), the situation is very different from 
the AV8' case. The channel coupling effect is rather modest, and the size of the 
1 Sq d+d elastic phase shift is already accounted for mostly in the d+d channel 
calculation. All these results are very consistent with the + spectrum in Fig. HI 

The large distortion effect of the deuteron clusters on the 1 S'o d+d scattering 
phase shift is expected to appear in the 3N+N phase shift as well because of 
the coupling between the 3N+N and 2N+2N channels. We display in Fig. [6] 
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Figure 6. So t+p elastic-scattering phase shift at energies below the d+d threshold. The 
solid line is a FULL calculation, while the dashed line is a 3N+N channel calculation. Crosses 
correspond to the i?-matrix analysis of Ref. [22] . 



the So t+p elastic-scattering phase shift at energies below the d+d threshold. 
The O2" state of 4 He is observed as a sharp resonance with a proton decay width 
of 0.5 MeV at about 0.4 MeV above the t+p threshold. The present energies 
(E r = 0.15 MeV for AV8', E r = 0.12 MeV for MN) calculated with a bound 
state approximation are slightly smaller than the experimental value, but they 
are consistent with a calculation (E r = 0.105 MeV and T/2 = 0.129 MeV for 
AV18+UIX, E r = 0.091 MeV and T/2 = 0.077 MeV for AV18+UIX+V^) with 
another realistic interaction (AV18) with three nucleon forces by Hofmann and 
Hale [22J. The calculated phase shifts appear slightly larger than that in the 
.R-matrix analysis (crosses in Fig. [6} [22]. It is noted that the phase shift changes 
so much even for a small change of the 0^" resonant pole position (~0.1 MeV) 
because it is very near to the threshold. The phase shifts in the FULL calculation, 
for both AV8' and MN cases, show a resonance pattern in a small energy interval 
and the overall energy dependencies of the phase shifts are similar to each other. 
However, the phase shifts obtained only in the 3N+N channel are quite different 
as indicated by the dotted lines in Fig. [6j In the case of the MN potential (right) 
the phase shift is already close to the FULL phase shift, while in the case of the 
AV8' potential (left) the phase shift is much smaller (by almost 90 degrees) and 
moreover shows no resonance pattern. 

By looking into the wave functions in more detail, we argue that the large 
distortion effect in the 1 So d+d and 3N+N coupled channels is really brought 
about by the tensor force. As shown in Table [2j the AV8' potential with TNF 
gives 5.8% and 8.4% (8.3%) Z?-state probability for d and t (h), respectively. 
Thus the d+d state in the 1 Sq state contains L = S = components (89%) as 
well as L = S = 2 components (11%), where L and S are the total orbital and 
spin angular momenta of the four-nucleon system. Similarly the 3N+N state 
in the 1 S'o state contains an L = S = component (92%) and an L = S = 2 
component (8%). Thus the tensor force couples both states with AL = 2 and 
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Figure 7. 1 Z?2 elastic-scattering phase shifts with the AV8' (left) and MN (right) potentials. 
The model space for the solid line is FULL. The t+p phase shift starts from the t+p zero energy 
and the d+d phase shift starts from the d+d threshold. The dashed line is the t+p phase shift 
with only the 3N+N channel and the dotted line is the d+d phase shift with only the 2N+2N 
channel. 



AS = 2 couplings, which are in fact very large compared to the central matrix 
element (AL = 0, AS* = 0). An analysis of this type was performed for some 
levels of 4 He in Refs. [39]. The MN potential contains no tensor force, so that 
the d+d and 3N+N channel coupling is modest. 

As listed in Table [H there are four channels, 5 52, 1 -t ) 2 5 3 ^2 and 5 D2, for 
J n = 2 + at energies around the d+d threshold. Among these states, we expect 
that the effect of the coupling between the 3N+N and 2N+2N channels occurs 
most strongly in 1 Z?2 as it appears in all physical channels. However, no sharp 
resonance is observed in 4 He up to 28 MeV of excitation energy, so that the 
coupling effect, if any, might be weaker than that observed in the 1 So case. 

Fig. [7] displays the 1 Z?2 elastic-scattering phase shifts obtained in three types 
of calculations, 3N+N (dashed line), 2N+2N (dotted line), and FULL (solid 
line). The t+p and d+d phase shifts start from the t+p (E c . m . = 0) and d+d 
thresholds, respectively. The phase shifts of the 3N+N and 2N+2N calculations 
are both slightly positive, indicating a weak attraction in the t+p and d+d inter- 
actions. In the FULL calculation, the t+p phase shift becomes more attractive 
and the d+d phase shift turns to be negative (repulsive). The present FULL 
calculation reproduces the calculation of Ref. [22] as expected. Though the effect 
of the coupling is slightly larger in the AV8' potential than in the MN potential, 
it is much less compared to the case of the 1 5o phase shift. This is understood 
as follows. In the 1 I?2 state, the main component of the wave function is given 
by the L = 2, S = state: Its probability is the same as that of 1 Sq, that is, 
92% in t+p and 89% in d+d. However, the probability of finding the state with 
L = 0, S = 2, which causes a strong tensor coupling, is more than one order of 
magnitude smaller than in the case of 1 Sq, namely 0.23% in t+p and 0.44% in 
d+d, respectively. The reason for this small percentage is that, to obtain L = 0, 
the incoming D-wave in the 1 L>2 channel must couple with the .D-components 
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in the clusters, but this coupling leads to several fragmented components with 
different L values. This relatively weaker coupling of the tensor force explains 
the phase shift behavior in Fig. [7J 
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Figure 8. Elastic scattering phase shifts of the 2 + state obtained in a FULL calculation with 
the AV8' (left) and MN (right) potentials. Solid line: 5 S 2 d+d phase shift, dotted line: 5 D 2 d+d 
phase shift, dashed line: 3 D2 t+p phase shift. 



In Fig. [8] we plot the t+p and d+d elastic-scattering phase shifts for other 
channels, 5 £2 (solid line), 3 I?2 (dashed line), and 5 L>2 (dotted line). We show 
only the FULL result, because the phase shifts with the truncated basis do not 
change visibly at the scale of the figure. The obtained phase shifts are not that 
different between the AV8' and MN potentials, and also consistent with the 
previous calculation [22]. Thus, the effect of the distortion of the clusters is very 
small for 2 + except for 1 D2- 

We have three channels for J n = 1 + , 5 D\, 3 D\ and 3 Si. No sharp 1 + res- 
onance of 4 He is observed experimentally up to 28 MeV of excitation energy. 
Another theoretical calculation neither predicts it [39], so that the coupling be- 
tween the 2N+2N and 3N+N channels is expected to be weak. Fig. [9] exhibits 
the t+p and d+d elastic-scattering phase shifts in the FULL calculation: 5 D\ 
d+d (solid line), 3 Z?i t+p (dashed line), and 3 Si t+p (dotted line). Only the 
FULL result is displayed because the phase shift change in other calculations is 
small. Both AV8' and MN potentials produce phase shifts quite similar to each 
other. 



5.3 Negative parity phase shifts 

As seen from Table HJ the main components of these negative parity states are 
considered to be 3 Pj. 

We compare in Fig. [10] the 3 P elastic-scattering phase shifts calculated with 
the AV8' (left) and MN (right) potentials. The truncated 3iV+iV (dashed line) 
and 2N+2N (dotted line) calculations are shown together with the FULL re- 
sult (solid line). The t+p phase shift of the 3N+N calculation is similar with 
both AV8' and MN potentials, while the d+d phase shift of the 2iV+2iV calcula- 
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Figure 9. Elastic scattering phase shifts of the 1 + state in a FULL calculation with the AV8' 
(left) and MN (right) potentials. Solid line: 5 Di d+d phase shift, dotted line: 3 Si t+p phase 
shift, dashed line: 3 Di t+p phase shift. 
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Figure 10. 3 Po elastic-scattering phase shifts calculated with the AV8' (left) and MN (right) 
potentials. See the caption of Fig. [T) 



tion behaves quite differently between the two potentials: the d+d phase shift is 
weakly attractive with AV8' but is very strongly attractive with MN. No typical 
resonance behavior shows up below the d+d threshold, which is in contradiction 
to experiment. In the FULL model that combines both 3iV+iV and 2N+2N 
configurations, however, the two potentials predict quite different phase shifts 
especially in the t+p channel. The t+p phase shift with AV8' becomes so attrac- 
tive that it crosses 7r/2, indicating a resonance at about 1 MeV above the t+p 
threshold. The d+d phase shift changes sign from attractive to repulsive. The 
result based on the AV8' potential is thus consistent with experiment. Further- 
more, we reproduce the flat structure of the 3 Po phase shift around several MeV 
above the t+p threshold which was discussed as the coupling to the h+n chan- 
nel [22]. On the other hand, the MN potential changes the t+p phase shift only 
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mildly and produces no sharp resonance behavior. The d+d phase shift changes 
drastically to the repulsive side. 

As seen in the above figure, the sharp 0~ state appears provided a full model 
space with a realistic potential is employed. The mechanism to produce this 
resonance is unambiguously attributed to the tensor force as discussed in Ref. [7] 
for the realistic interaction G3RS [41j . According to it, the 0~ state consists of 
only two components, L = S = 1 (95.5%) and L = S = 2 (4.5%), ignoring a tiny 
component with L = S = 0. The L = S = 2 component arises from the coupling 
of the incoming P-wave with the -D-states contained in the 3N and d clusters. 
All the pieces of the Hamiltonian but the tensor force have no coupling matrix 
element between the two components. The uncoupled Hamiltonian thus gives a 
too high energy to accommodate a resonance. The tensor force, however, couples 
the two components very strongly, bringing down its energy to a right position. 

The second lowest negative parity state has spin-parity 2~. The physical 
channel for this state is only 3 i-2 as seen in Table [TJ Fig. [11] compares the 3 P2 
elastic-scattering phase shifts in a manner similar to Fig. [TOl The phase shift 
obtained with the MN potential is almost the same as the 3 Po phase shift, which 
is consistent with the previous result [39] that the energies of the negative parity 
states calculated with the MN potential are found to be degenerate. In the case 
of the AV8' potential, the 3 P2 phase shifts grows significantly in the FULL cal- 
culation, indicating a resonant behavior. The coupling effect between the 3N+N 
and 2N+2N channels is however much less compared to the 0~ state. This is 
because the incoming P-wave coupled to the -D-states in the clusters gives rise 
to several L values to produce the 2~ state and therefore the tensor coupling 
does not concentrate sufficiently to produce a sharp resonance. 
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Fig. [12]displays the 3 Pi and 1 P\ elastic-scattering phase shifts calculated with 
the AV8' (left) and MN (right) potentials. Note that no physical d+d channel 
exists in the case of the 1 P\ state. Because both FULL and 3N+N calculations 
give almost the same phase shifts, only the FULL result is shown in the figure. 
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The 3 .Pi phase shift calculated with the MN potential is again almost the same as 
those of the 3 Pq and 3 P2 cases, supporting that the three negative parity states 
become almost degenerate. The 3 P\ elastic-scattering phase shift calculated with 
the AV8' potential is qualitatively similar to that of 3 i-2- The attractive nature of 
the t+p phase shift becomes further weaker, and to identify a resonance appears 
to be very hard. Even though it is possible in some way, its width would be a 
few MeV, which is not in contradiction to experiment. The l P\ phase shifts are 
very small in both AV8' and MN cases. 

For the negative parity states, the FULL model with the AV8' potential gives 
results that are consistent with both experiment and the theoretical calculation 
of Ref. [39] • We have pointed out that the phase shift behavior reveals the impor- 
tance of the tensor force particularly in the case of 0~ . Its effect is often masked 
however by the coupling between the D states in the clusters and the incoming 
partial wave. 
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In this subsection, we investigate the phase shifts of the negative parity states 
which have dominant T = components. In Fig. [T3l we represent three exper- 
imental negative parity T = energies (left). The states are observed at —7.29 
(0~), -6.46 (2") and -4.05 (1") MeV below the four-nucleon threshold @Dj. The 
former two are located below the d+d threshold and their widths are 0.84 and 
2.01 MeV, respectively, whereas the last one is above the d+d threshold and its 
width is fairly broad (6.1 MeV). Here, we calculate these energies as —7.57 (0~), 
—6.82 (2~) and —5.95 (1~) MeV, which are approximated by the half-value po- 
sition from the maximum phase shifts. The present calculation is not projected 
out to T = 0, but the dominant configurations of t+p, h+n and d+d elastic 
scattering are T = 0. Our calculated energies with AV8' reproduce the ordering 
of _ , 2^ and 1~ (middle in Fig. 1 13[> . The splitting between the two lower states 
_ and 2— is reproduced, but the experimental 1~ energy is higher than the 
calculation. However, the determination of the energy for such a high energy 
state with large decay width (6.1 MeV) is very difficult from both experimental 
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and theoretical side, and it usually has a large ambiguity. 

This type of analysis was done by Horiuchi and Suzuki, who applied the 
correlated Gaussian basis with two global vectors to study the energy spectrum 
of 4 He [39] . Because the results of Ref . [39] are based on approximate solutions 
that impose no proper resonance boundary condition, it is interesting to see how 
the tensor force changes the phase shifts in the negative parity states as shown 
in this subsection. These authors also found that the negative parity states with 
T = turn out to be almost degenerate when the MN potential that contains 
no tensor force is employed. In the present calculation, three states (0~, 1 _ , 2 _ ) 
completely degenerate at the same energy, E = —6.64 MeV (right), and the same 
phase shift pattern (solid lines in Figs, [TO l [TT | [T2l . Thus we can expect to see a 
clear evidence for the tensor force in the scattering involving the negative parity 
states. 



6 Summary and conclusion 

We have investigated the distortion of clusters appearing in the low-energy d+d 
and t+p elastic scattering using a microscopic cluster model with the triple 
global vector method. We showed that the tensor interaction changes the phase 
shifts very much by comparing a realistic interaction and an effective interaction. 
In the present ab-initio type cluster model, the description of the cluster wave 
functions is extended from a simple (0s) harmonic-oscillator shell model to a 
few-body model. To compare distortion effects of the clusters with realistic and 
effective interactions, we employed the AV8' potential as a realistic interaction 
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and the MN potential as an effective interaction. 

For the realistic interaction, the calculated 1 5o phase shift shows that the 
t+p and h+n channels strongly couple with the d+d channel. These channels 
are coupled because of the tensor interaction. On the contrary, the coupling of 
these 3N+N channels plays a relatively minor role for the case of the effective 
interaction because of the absence of tensor term. In other words, the 3N+N 
channels strongly affect the d+d elastic phase shift with the realistic interaction, 
but not with the effective interaction. 

For the 2 + phase shifts, there is a 1 D 2 component in all physical channels 
(d+d, t+p and h+n). The coupling of the 2N+2N and 3N+N channels in 5 D 2 
is weaker than in 1 S 2 because of a weaker tensor coupling as discussed in section 
El and the calculated phase shifts are very similar for the realistic and effective 
potentials. For other positive parity cases, the phase shift behavior of the realistic 
and effective potentials are very similar, and the coupling between the 2N+2N 
and 3N+N channels can be neglected or is very small. Furthermore, the tensor 
interaction makes the energy splitting of the 0~, 2~ and 1~ negative parity states 
of 4 He consistent with experiment. No such splitting is however reproduced with 
the effective interaction. 

We believe that the physical picture obtained in the large model space with 
the realistic interaction should be close to the real physical situation. It is needless 
to say that ab-initio reaction calculations are very important to understand the 
underlying reaction dynamics involving continuum states. Simpler calculations 
using effective interactions in the same framework, as carried out in the present 
paper, are also meaningful because we can understand more clearly the effect 
of the tensor force by comparing both calculations. The reaction calculations 
with the microscopic cluster model, whose model space and interactions are 
restricted, have been successfully applied to many heavier nuclei. Therefore, it 
is instructive to see the difference from the realistic interaction by employing a 
simple conventional effective interaction as MN in the few-nucleon systems. 

It will be quite interesting to see the importance of the tensor force in 
reaction observables of four nucleons. As a direct application of the present 
study the radiative capture reaction d(d, 7) 4 He at energies of astrophysical 
interest is of prime importance. It is expected to take place predominantly via 
E2 transitions [121 ETJ 031 [201 [9] . As is seen from Table [U the two deuterons 
can approach each other in the S- wave only when J w is either + ( 1 S'o) or 
2 + ( $2)- The former case is excluded because a radiative capture reaction of 
+ — > + is forbidden in the lowest-order electromagnetic interaction, and hence 
the E2 transition should be predominant. If there were no tensor force present, 
the radiative capture would be suppressed near E = because neither d nor 
4 He would have a D-w&ve component in contradiction with the flat behavior 
of the astrophysical S- factor |44| . The tensor force strongly changes this story 
because it can couple S- and D-waves, bringing a significant amount of D-state 
probability in both 4 He and d. Details of this analysis will be reported elsewhere. 
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Appendix A: Definitions of recoupling coefficients 

We define an auxiliary coefficient Z that appears in the coupling 

[[Ka(e 1 )[Yi,(ei)y6(e2)]o]a -► Z{ab)[Y a+b { ei )Y b {e 2 )] a . 
By introducing a coefficient 



C(ab, c) = 



(2o + l)(26 + l) 



a b 0|c 0} 



4tt(2c + 1) 

for the coupling [Ya(ei)Yf,(ei)] c = C(ab,c)Y c (ei), we can express Z as 

1 



Z(ab) = 



2(a + b) + 1 
(2o + l)(26 + l) 



C(ab,a + b) = —={a b 0\a + b 0) 



47T 



Note that C(ab, c) vanishes unless a + b + c is even. 

The coefficients that appear in Sect. [4] are given as follows: 

X(a b c) = Z(ab)Z(ac)C(bc, b + c)U(a+c c a + b b; a b+c) 

R 3 {a bc) = Z(ab)Z(a + b c), 

W(a bcqQ, a b' c q Q' , £ t) 

q c Q a b q 

q c Q a b q 

I' c+c' I a + a' b+b' £' 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
(A.5) 



C(aa ,a+a')C(bb' ,b+b')C{cc , c+c). (A.6) 



Appendix B: Matrix elements for various operators 

The purpose of this appendix is to collect formulas for various matrix elements. The main 
procedure to derive the formulas is sketched in Sect. 0] More details for the case of two global 
vectors are given in Ref. [?]• 



The total kinetic energy 



B.l Kinetic energy 

Let 7Tj denote the momentum operator conjugate to Xj, ttj = — ih- 
operator for the A^-nucleon system with its center of mass kinetic energy being subtracted takes 
the form 

N - 2 -a i 

nN - (B.l) 



— — 7rA7r 

2m 2A^m 2 



where 7Tjv = 5Z^=i Pi i s the total momentum, tt = (7ri,7T2, . . . , 7Tjv-i), and A is an (N — 1) x 
(N — 1) symmetric mass matrix. Defining N — 1-dimensional column vectors I\ as 



Y i = A'B- 1 Ui (i = l,2,3), 
r i = -AB- 1 u i (i = 4,5,6) 



(B.2) 
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and an (TV — 1) x (N — 1) matrix Q 

Q lJ =2f l Ar j , (B.3) 
we can calculate the matrix element for the kinetic energy through the overlap matrix element 

<F'\l*An\F) = ^ (V|> f ^)<F'|F>, (B.4) 

where 

R = 3Tr{B~ 1 A'AA). (B.5) 
The pij values are defined in Eq. (|4.11[) . 

B.2 5- function 

A two-body interaction V(ri — rj) can be expressed as 

V(n -rj) = j drV{r) 5(n - rj - r). (B.6) 

Once the matrix element of S(ri — rj — r) is obtained, the matrix element of the interaction is 
calculated by integrating over r the ^-function matrix element weighted with the form factor 
V(r). Similarly, for a one-body operator 

D(ri — x N ) = j drD(r) 5(r 4 - x N - r), (B.7) 

its matrix element can be obtained from that of the S- function. Because both ri — rj and ri — XN 
can be expressed in terms of a linear combination of the relative coordinate Xi, it is enough 
to calculate the matrix element of S(wx — r), where w = (wi,W2, . . . , wn-i) is a combination 
constant to express Vi — r j or r, — xn ■ 

The matrix element of the 5-function is given by 

(F'\8(wx - r)\F) 

■(2n) N ' 1 \ 3/2 (Jj \ (-l) i i+^+^v / 2TTT / c \S/a 



n J V2.L + 1 /_c_ 

detB J VfJ L ' ) V2i' + 1 V2tt/ 

j2(lm^\l'm')y:,(?) y: ( ri(-^r 



f_iY e +* 

x 22 - I U(LLk£';£L') Ojpj-Jwhstt'K) 

i 12 i i5 u'T 12 L i5 T y (2£ + 1)(2L + 1) 

x W{p 1 p 2 pzii2Z,L 1 -p 1 L 2 -p 2 L 3 -p 3 L 12 L,LL 12 ) 
x Wip^PhPd^ , Li-pi L 5 ~p 5 L 6 ~p 6 L45L, L'Las) 



™ii \i<3 



0(mj\ Li— pi L 2 ~P2 Lz—pz £4— £>4, £5— Ps L & —p e , L^LasL), (B.i 



with 



c=(wB 1 w) \ 7i=ioB 1 u i , Pij = pij-cjijj. (B.9) 
The summation over non-negative integers ny and pi is restricted by the following equations 

"12 + "13 + "14 + "-15 + ma + pi = Li, 

ni2 + "23 + "24 + "25 + "26 + P2 = L 2 , 
"13 + "23 + "34 + "35 + "36 + PZ = 1/3, 
7114 + "24 + "34 + "45 + "46 + Pi = La, 
"15 + "25 + "35 + "45 + "56 + Pb = £-5, 

"16 + "26 + "36 + "46 + "56 + P6 = ^6- (B.10) 



30 



Here 0(pi; IviIa^U! k) is defined as a coefficient that appears in the coupling of a product of six 
terms 

6 

ll[Y Pi (ei)Y Pi (r)] 00 = ]T 0(p i; 4^45^' '«) 

x [[[[ii I (ci)y„(e a )]/ I!i yp3(es)] < [[y M (e4)y p5 (e5)]^ 5 K P6 (e 6 )]^] K y K (f)]oo, (B.ll) 
and it is given by 



C(pip 2 , £l 2 )C*( W^)C*(P4P5^45)C , (£ 4 5P6^')C(«', «)• (B.12) 



nLi(2^ + 1) 

The r-dependence of the matrix element (|B,8|I is 

r Pl+ P 2+P3+P4+P5+P6y- K . (B.13) 



For a central interaction, V(r) is a scalar function, and the sum over k in Eq. ()B.8[) is limited 
to 0. For a tensor interaction, the angular dependence of V(r) is proportional to Y2(r), and k 
is limited to 2. The electric multipole operator is a special case of one-body operator, so that 
one can make use of the formula (|B.8|) to calculate its matrix element. More explicitly, we give 
the matrix element of V {\wx\)Y Kfl {wx) that includes all the cases mentioned above: 



F'\V{\wx\)Y^{wx)\F) 

(2tt) jv - 1 \ 3/2 / A „ \ (-l) il+i2+L 3 v / 2TTT 



detB J \±± 'J V2L' + 1 



+P2+P3+P4+P5+P6 



(_1)W 



x = U{LLk£';£L') Q( Pi ;l 12 £ i5 ££'K) 

i^t^u'L^L^L \] {21 + 1)(2L + 1) 

x W(pip2P3£i2£, Li— pi L2—P2 L3—P3 L\2h,LL\2) 
x W (pAPsP^Abt' , L A -p A L 5 -p 5 L 6 ~p 6 LasL,L'Las) 



x 0{mj\ L\— pi Li— pi L 3 —p 3 La— p4,L 5 —p 5 La — P6,Li 2 LasL), (B.14) 
with the integral of the potential form factor 



(c) = (^) c"y drr n+m V{r)e-^ cr2 . (B.15) 

In case V{r) takes the form of r q e~ pr2 ~ p ' r {q > -m), the integral l ( n m) (c) can be obtained 
analytically, giving a closed form for the matrix element. 

It should be noted that the matrix element for a special class of a three-body force can be 
evaluated with ease. For example, if the radial part of the three-body force has a form 

Vtnf = exp(-pi(r, ; - r^f - p 2 {rj - r k )' 2 - p 3 (r k - r^) 2 ), (B.16) 

the exponent can be rewritten as — x£lx with an (TV — 1) x (N — 1) symmetric matrix O = 
piWijVHj + p2W jk wJ k + p3W ki Wh~i, where Wij, w jk and w ki are defined by n — Vj = wtjX, 
fj — r k — wJkX and r k — ri — WkiX. Thus the matrix element reduces to that of the overlap 
with A being replaced with A + 20 

(F'\Vtnf\F) = (F'\F LlL2(Ll2)L3LM {ui,u 2 ,u 3 ,A + 2n,x)}. (B.17) 
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B.3 Spin-orbit potential 

The spatial form of a spin-orbit interaction reads 

VQn-TiWn-ri) x (B.18) 

where (a x b) M = 0, ±1) stands for the pith component of a vector product of a and b. As in 
the 5-function matrix element, the spin-orbit potential is written as 

V(\wx\)(wx x Ctt)m. (B.19) 

where \{p i — Pj) is expressed in terms of the momentum operators n, <^7r = Y^i^ 1 O^i- 
The matrix element of the spin-orbit potential is given by 



F'\V(\wx\)(wx x (-k)„\F 

4ny/2h f(2n) N - 1 \ 3/2 f TT Bl \ (-l) il+i2+i V 



detB J yi = l 'J V2L' + 1 

{LMML'M')^ (fli-^r^^j ^ +P2+P3+Pi+P5+Pe (c) 



Pi \i=l 

\l+7 



x Y, J 1 _ =U(LLU;£L') 0{p i ;l^U^(!\) 

li_il ih ll r lizl i {i I'l 12 l 45 l y (2^ + 1)(2L + 1) 

x ^(cr) fc r fc ( Pl ,£ 12 f45^'J 12 ^ 5 II')E (fl(-p y )" y ^^) 

k — l n»j \i<j ™ l 3 / 

x 0{nif,L\—px L 2 — p 2 L 3 — p 3 La — Pa L 5 —p 5 L 6 — pq,L 12 LasL), (B.20) 
where p\ (k — 1, 2, . . . , 6) is 

=Pi + (B.21) 
and where the non-negative integers riij and pi are constrained to satisfy the equations 

nu + ni3 + tiia + nis + nie + Pi = L lt 
ni2 + n 23 + n 2 A + n 25 + n 26 + p 2 = L 2 , 
ni3 + n 23 + n 3 A + n 35 + n 36 + p 3 = L 3 , 

TilA + TL2A + n-3A + «4S + ™46 + Pa = La, 

nis + n 25 + n 35 + n 45 + n 56 + p\ = L 5 , 

nia + n 26 + n 3(i + n 46 + n 56 + pg = L 6 . (B.22) 
The symbol (Cr)fc stands for the factor 

6 

(Cr)* = Z)C<(r fc )<> (B.23) 



where (Tk)i is the ith element of the column vector Tk defined in Eq. (|B.2[I . The coefficient Tk 
appears in the coupling 

[Yx(e k ) [[[Y P1 (ei)Y M {e 2 )] ll2 Y P3 (e 3 )] t [\Y P4 (e 4 )Y P5 (e 5 )] ti5 Y P6 (e 6 )Mi]i„ 

x [[[Y p * (er)^ (e 2 )] Il2 Y p * (e 3 )] ? [[Y p * (e 4 )Y p g (e 5 )b 45 Y p| (e«)]y]i„. (B.24) 
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The coefficients Tk(pi, £i 2 lis££' , £i 2 £a$£ £ ) are given below: 

Ti = U{l£U'f£l)U{U i2 £p 3 -li2£)U{lpil i2 p2-,Px+l £ia)C(lp l5 pi+l) 

r 2 = -(-i)' 12+ " 12 c/(m/;n)[/(w 12 Ip3;^)i7(Wi2Pi;p2+i ^ 2 )C(ip 2 ;p2+i) 

T 3 = -(-l) i+I U(im';£l)U(lp s l£ir,P3+U)C(lp 3 ;p3+l) 

Ti = (-1) £ ' +? [7(1/W;/I)f/(W45I'P4;I45/)(7(1P4?45P5;P4 + 1 ^4 5 )C(lp4; P4 + 1) 

T 5 = -(-l) f ' +?+ ^ +745 C/(l^l^Z'l)[/(W45l'p6;l4 5 /)C/(lp 5 ^4 5 P4;P5 + l 4 5 )C(lp 5 ; P5 + 1) 
T 6 = -(7(l/W;/l)?7(lp 6 ^45;P6 + l /)C(lp 6 ;p 6 + l). (B.25) 

Appendix C: Momentum representation of correlated Gaussian basis 



The Fourier transform of the correlated Gaussian function (|3.9|) defines the corresponding basis 
function in momentum space. The momentum space function is useful to evaluate those matrix 
elements which depend on the momentum operator [7j. Suppose that we want to evaluate the 
matrix element of a two-body operator V(p i — p-) or a one-body operator D(p i — -^ttn)- 
Obviously evaluating the matrix element can be done more easily in momentum space. For 
this purpose we need to obtain the Fourier transform of the coordinate space function. A 
great advantage in the correlated Gaussian function F is that its Fourier transform is a linear 
combination of the correlated Gaussian functions in the momentum space. Thus by expressing 
p i — Pj or p i — jfTTN as <^7r, we can calculate the matrix element of the momentum-dependent 
operators in exactly the same way as in the coordinate space. 

As in the case with two global vectors [7], the transformation from the coordinate to 
momentum space is achieved by a function 

$(fc,a;) = exp(ifca;), (C.l) 

(27T)i( Ar - 1 ' 

where fc is an (N— l)-dimensional column vector whose ith element is fci. With a straightforward 
integration together with the recoupling of angular momenta, we can show that 

(<6(fe, x)\F LiL2 ( Li2 ) L3 lm(ui,U2,U3, A, x)) 

(_j\L 1 +L 2 +L 3 



(detAW 2 E 'C(LiL 2 (L 12 )L 3 L;e 1 e 2 e 3 £ 12 ) 

X-Pii^i^c 2 hri-L-t* (iii) L 3 -e 2 -e 3 LM{A~ 1 ui,A~ 1 u 2 ,A~ 1 u 3 ,A~ 1 ,k), (C.2) 
where the coefficient K, is given by 
K,(LiL 2 (Li 2 )L 3 L;£ie 2 £ 3 £i 2 ) 

_ (^L-Ls+ia+ts-tii B Ll B L2 B L3 



y / 2L+T Be 1 Be 2 Bi 3 Bi J1 _e 1 -e 2 Bl 2 -i 1 -i 3 BL 3 -e 2 -i 3 

x y /(2£ 1 + l){2£ 2 + l){2£ 3 + l){2(L 1 ^£ 1 -£ 2 ) + l)(2(L 2 ^£ 1 -£ 3 ) + l)(2(L 3 -£ 2 -£ 3 ) + 1) 
xX{£ 1 £ 2 £ 3 )Z(L 1 -£ 1 -£ 2 £i+£ 2 )Z(L 2 -£ 1 -£ 3 £ 1+ £ 3 )Z(L 3 -£ 2 -£ 3 £ 2 +£ 3 ) 

" L x -£i-£ 2 Li £-,+£ 2 " 



x I7(^i2 L 3 -£ 2 -£ 3 L12 L 3 -L£ 2 +£ 3 ) 



L 2 -£t-£ 3 L 2 h+£ 3 
£12 L 12 £ 2 +£ 3 



x^A-^Y 1 (mA -1 ^)* 2 (uiA^uaf 3 , (C.3) 

where Z and X are defined in Appendix A. Non-negative integers £i run over all possible values 
that satisfy £1+^2 < L\, £i+£ 3 < L 2 , £ 2 +£ 3 < L 3 . The value of £± 2 is restricted by the triangular 
relations among (£ i2 , Li — £1 — £ 2 , L 2 — £\ — i 3 ) and {£\2,Lr 2 ,£ 2 +£ 3 ). 
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